Non-adiabatic dynamics of two strongly coupled nanomechanical resonator modes 



Thomas Faust, Johannes Rieger, Maximilian J. Seitner, Peter Krenn, Jorg P. KotthausQand Eva M. Weig[[ 

Center for Nano Science ( CeNS) and Fakultdt fur Physik, 
Ludwig-Maximilians-Universitat, Geschwister-Scholl-Platz 1, Miinchen 80539, Germany 

The Landau- Zener transition is a fundamental concept for dynamical quantum systems and has 
been studied in numerous fields of physics. Here we present a classical mechanical model system 
exhibiting analogous behaviour using two inversely tuneable, strongly coupled modes of the same 
nanomechanical beam resonator. In the adiabatic limit, the anticrossing between the two modes 
is observed and the coupling strength extracted. Sweeping an initialized mode across the coupling 
region allows mapping of the progression from diabatic to adiabatic transitions as a function of the 
sweep rate. 
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The time dynamics of two strongly coupled harmonic 
oscillators follows the Landau-Zener model [T]-[4] , which 
is used to describe the quantum mechanical mode tun- 
neling in a non-adiabatic transition. This phenomenon 
is observed and utilized in many areas of physics, e.g. 
atomic resonances [5] , quantum dots [6 , superconducting 
qubits [7 1 and nitrogen- vacancy centers in diamond [8] . It 
is also possible to create classical model systems exhibit- 
ing the same time evolution, which until now have been 
restricted to optical configurations [9j [10] . Such systems 
are well suited for the study of diabatic behaviour over 
a wide parameter space; for example nonlinearities could 
be readily introduced, potentially leading to chaotic be- 
haviour [11] • 

Nanomechanical resonators with frequencies in the 
MHz range can be realized with high mechanical quality 
factors [I2j [13] and easily tuned [14] in frequency. This 
makes them particularly well-suited for exploration of 
their coupling to other mechanical, optical or electrical 
microwave resonators. Strong cavity coupling in the op- 
tical or microwave regime has been widely studied as it 
enables both cooling and self-oscillation of the mechan- 
ical modes [T5HT81 . In addition, the time-resolved Rabi 
oscillations between a strongly coupled two-level system 
and a micromechanical resonator have been observed [19 . 

Purely mechanical, static coupling between different 
resonators [20 -23 and between different harmonic modes 
of the same resonator [24 has also been demonstrated. 
Here, we explore the coupling between the two fundamen- 
tal flexural modes [25 of a single nanomechanical beam 
vibrating in plane and out of plane, respectively. We 
study the adiabatic to non-adiabatic transitions between 
the two strongly coupled classical mechanical modes in 
time-dependent experiments, in correspondence to the 
Landau-Zener transition. 

The nanomechanical high stress silicon nitride string 
used in this work is shown in Fig. [T] Two parallel gold 
electrodes vertically offset to the beam are used to di- 
electrically couple the beam oscillation to an external 
microwave cavity with a quality factor of « 70 at a res- 
onance frequency of 3.44 GHz [26 . Displacement of the 




FIG. 1. (color online). The SEM micrograph of the 55 um 
long and 260 nm wide silicon nitride string resonator (green) 
taken at an angle of 85° also depicts the two adjacent gold 
electrodes (yellow) used to dielectrically drive, tune and read 
out the resonator motion. The arrows denote the two me- 
chanical modes, one oscillating parallel and the other perpen- 
dicular to the plane of the chip. The simplified measurement 
scheme [26] shows the connection of the electrodes to the read- 
out cavity (gray box) and the microwave bypass capacitor in 
the bottom left. 



resonator leads to a change in capacitance between the 
two electrodes, thereby detuning the resonance frequency 
of the microwave circuit and creating sidebands with a 
frequency offset equal to the mechanical eigenfrequency. 
The inductively coupled microwave cavity is driven by a 
signal generator; the transmission signal is demodulated 
and fed to a spectrum analyzer as depicted in Fig. [I] and 
described in more detail in [26 . In addition, a microwave 
bypass capacitor is used in the ground connection of one 
electrode which allows application of additional DC bias 
and RF voltages to the electrodes. This is used to ac- 
tuate the mechanical resonator via the dielectric driving 
mechanism [14 , 27]. At the same time, the dielectric cou- 
pling provides a way to tune the resonance frequency 
of the two mechanical modes: The static electric field 
between the electrodes polarizes the dielectric resonator 
material which is then attracted to high electric fields, 
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FIG. 2. (color online). Both mechanical modes can be tuned 
in opposite direction by increasing the DC bias voltage applied 
to the electrodes. The signal power of the driven resonances 
is shown color-coded versus DC voltage and drive frequency. 
Note the clear avoided crossing between the two modes. The 
three circles denote the initial state (I) and two possible fi- 
nal states after an adiabatic (A) or diabatic (D) transition 
through the coupling region, as described in the text. 



thereby changing the spring constant of the modes via 
the resulting force gradient [14] . In the chosen geome- 
try, where the bottom of the electrodes is flushed with 
the top of the beam [26], a rising DC bias voltage causes 
the frequency of the in-plane mode to decrease and the 
out-of-plane frequency to increase [28] . All experiments 
are performed at room temperature at pressures below 
5 x 10 -4 mbar. 

At low DC bias voltages, the in-plane mode of the 
55 um long beam has a higher resonance frequency than 
the out-of-plane mode. This is a result of the 260 nm 
beam width exceeding the beam's thickness of 100 nm, 
which leads to a higher rigidity for the in-plane mode [13] . 
Thus, by increasing the DC bias voltage, we are able to 
tune the two modes into resonance at a common fre- 
quency of approximately 6.63 MHz. The coupling be- 
tween the modes has been observed for several resonators 
on various chips and is at least partially caused by the 
spatially inhomogenous electric field [29] . There might 
also be an additional, purely mechanical coupling medi- 
ated by the prestress in the beam. The characteristic 
avoided crossing diagram of two coupled oscillators can 
be obtained by measuring the driven response of the two 
modes at different DC bias voltages, as shown in Fig. [2] 

Splitting this diagram into an upper and lower branch 
and fitting each dataset with a Lorentzian allows the ex- 
traction of the resonance frequencies and quality factors 
for each DC bias voltage applied to the electrodes. Both 
modes exhibit a quality factor of approximately 80 000, 
somewhat lower than in previous measurements [26 , pre- 
sumably caused by fabrication imperfections. The eigen- 
frequencies extracted from the anticrossing diagram are 



_6 

N 

c 6 




66 
65 
64 
63 
62 
61 
60 
6.59 



IN ™ 







■ 


T7?ir 


■ 
■ 


: out 




IN ■ 


- i 






4 


5 6 7 8 
DC voltage (V) 


9 



FIG. 3. (color online). Frequency of the upper (red) and lower 
(blue) branch versus DC bias voltage. Each dot represents 
a value extracted from a Lorentzian fit of the data shown 
in Fig. [2] the solid lines are a fit of the theoretical model 
described in the text. IN and OUT denote the in- and out- 
of-plane mode of the beam. 

depicted in Fig. [3] A few data points around 6.5 and 
7.4 V in the upper branch were omitted because of an 
insufficient signal to noise ratio. 

For our system, the standard model of two coupled 
harmonic oscillators [30] needs to be expanded, as both 
oscillators react differently to the tuning parameter (the 
DC bias voltage). We use the generalized differential 
equation for the displacement u n of each mode n (n = 
1,2) 



with 



k n 



k q — \~ k(2 



(i) 



(2) 



where m e fr denotes the effective mass and 7 = uo/Q the 
damping constant of the resonator (identical for both 
modes), k c the coupling between the two modes and k n 
the spring constant of the respective mode. As the DC 
bias voltage polarizes the resonator material and creates 
an electric field gradient, the additional force gradient 
seen by the beam depends on the square of the voltage. 
We use a second-order series expansion around Uo to de- 
scribe the tuning behaviour: k n — ko + K n (U — Uq) + 
^n(U — Uq) 2 with n n and A n as linear and quadratic tun- 
ing constants, assuming that both modes have the same 
spring constant k$ at the voltage Uq corresponding to 
zero detuning. Note that the influence of the quadratic 
term is less than 15 % in the whole voltage range [29 . The 
two solutions of the differential equation describe the 
two branches, their fit to the experimental data is shown 
as solid lines in Fig. [3] The extracted frequency splitting 




7.77 kHz 
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at £/o=6.547V is much larger than the linewidth of 
7/27T = 83 Hz, thus the system is clearly in the strong- 
coupling regime. 

When slowly (adiabatically) tuning the system 
through the coupling region, the system energy will re- 
main in the branch in which it was initialized, thereby 
transforming an out-of-plane oscillation to an in-plane 
motion (and vice- versa for the other mode). At high tun- 
ing speeds, the diabatic behavior dominates and there 
is no mixing between the modes. This classical be- 
haviour [30] is analogous to the well-known quantum me- 
chanical Landau-Zener transition. The transition prob- 
abilities are identical in the quantum and classical case: 
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where the change of the frequency difference between the 
two modes in time 



d(oui - 002) 

a = — using u n 

ot 



kp + k n 



(5) 



denotes the tuning speed [29] . 

The measurement sequence is depicted in Fig. |4ja): 
The system is initialized at point I (see Fig. [2| by apply- 
ing a 6.6647 MHz tone and a DC bias voltage of 3.6 V to 
the electrodes. At t = 0, the voltage (blue line) is now 
ramped up to 9.1 V within time r. As the start and stop 
frequencies are kept constant throughout the experiment, 
changing r changes the tuning speed a and therefore the 
transition probability. Thus, the system's energy is dis- 
tributed between point A or D (see Fig.[2|, depending on 
the ramp time r. At t = 0, the mechanical resonator gets 
detuned from the constant drive frequency. Therefore its 
energy starts to decay as reflected by the decreasing sig- 
nal power (green dashed line in Fig. [4ja)). After a short 
additional delay of 5 (to avoid transient artifacts in the 
measurement), the decay of the mechanical oscillation is 
recorded with the spectrum analyzer. An exponential fit 
to the signal power, symbolized by the dotted black line 
in Fig. |4ja) allows the extraction of the oscillation mag- 
nitude at t = r which is normalized to the magnitude 
measured before the transition at point I to account for 
slight variations in the initialization. This experiment is 
repeated with many different ramp times r and with the 
detection frequency of the spectrum analyzer set to mon- 
itor either point A or D. The results of these measure- 
ments are shown in Fig.Qb). The data clearly shows the 
expected behaviour: For short ramp times below 0.2 ms, 
the diabatic behaviour dominates. For long ramp times, 
the adiabatic transition prevails, even though mechanical 
damping decreases the signal for large r. 

As can be seen in the inset of Fig. [4^b), the sum of 
the two curves perfectly follows the exponential decay of 
the mechanical energy (solid line). This decay in ampli- 
tude between t = and t = r has to be accounted for in 
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FIG. 4. (color online). The measurement sequence of the 
time-resolved experiment is shown in (a). At t = 0, the DC 
bias voltage (blue line) is ramped up in the timespan r, after 
the delay 5 the measurement of the mechanical signal power 
(green dashed line) starts at point A or D in Fig. [2] and a fit 
(black dotted line) is used to extract the magnitude of the 
beam oscillation at t = r. The normalized signal power at 
t — r and thus the transition probability obtained for differ- 
ent ramp times r measured at point A in Fig. [3] (blue trian- 
gles) or point D (red dots) is plotted in (b) together with the 
theoretical model described in the text (solid lines). The in- 
set shows the sum of both measurements and displays a clear 
exponential decay. The corresponding decay probability is 
represented by a green dashed line in the main plot. 



the theoretical model and therefore an additional decay 
term e~ 7t is introduced to equations Q . The solid lines 
in Fig. [4^b) show the resulting transition probabilities to 
point A and D and are calculated by using the a and T 
obtained from the data in Fig. [3] The measured data was 
rescaled by a constant factor with no free parameters to 
represent the probability distribution of the resonator's 
energy after a transition [29] . A third state, represent- 
ing the probability that the mechanical energy decays, 
is required to keep the sum of the probabilites at one. 
It is determined from the inset and shown as a dashed 
green line in Fig. [4j The corresponding decay constant 
I/7 = 1.92 ms is identical to the one extracted from the 
spectrally measured quality factor. Note that dynamics 
with a time constant much smaller than I/7 are observed, 
demonstrating coherent control of the system. 

In conclusion, we utilize the strong coupling between 
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two orthogonal modes of the same nanomechanical res- 
onator by tuning these two modes into resonance to an- 
alyze their time-dependent dynamics. After characteriz- 
ing the coupling, we are able to model the time-resolved 
transition behaviour between the two modes. The entire 
dynamic range between fast and coherent diabatic and 
slow adiabatic passages is accessible in the experiment, 
and a good agreement between theory and experiment is 
observed. 

The experiment is conducted with approximately 
10 9 phonons in the vibrational mode of the resonator, 
thus not the single-particle probability function but the 
energy distribution of the ensemble is measured in the 
classical limit. Since the (strongly) nonlinear regime 
of the utilized nanomechanical resonator can be easily 
accessed, the presented system could also be used to 
study the coupling and the time-dependent transitions of 
two nonlinear oscillators [3TJ [32] and the development of 
chaotic behaviour [11] [32] in the classical regime. Com- 
bining cavity-pumped self-oscillation [16] with the cou- 
pled resonator modes presented here allows the study of 
synchronisation and collective dynamics in nanomechan- 
ical systems, as theoretically predicted [33), [34] . Further- 
more, after the recent breakthrough in the ground state 
cooling of mechanical resonators [17-19 , the coupling be- 
tween two quantum mechanical elements becomes acces- 
sible. 
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SUPPLEMENTAL MATERIAL TO " NON- ADIAB ATIC DYNAMICS OF TWO STRONGLY COUPLED 

NANOMECHANICAL RESONATOR MODES" 



Modeling and data analysis 

Modeling the anticrossing 
The two coupled differential equations for the displacements u n (n=l,2) are 



m e ff-^- +m e ff7i-^ + hui + k c (ui -u 2 ) = (SI) 

d 2 u 2 du2 i / \ /r , x 

m e ff-^- + m e ff72-^ + fc2^2 + fc c (^2 - ^i) = 0, (S2) 

using the effective mass of the resonator (which, for a doubly clamped beam, is half its total mass) m e ff, the spring 
constants fc n , the coupling constant fc c and the two different damping constants j n . They are solved using an oscillatory 
ansatz u n = a n e lUJt . This leads to 

— cj 2 ra e ff<2i + fciai + k c (ai — a 2 ) + im e ffCJ7iai = (S3) 
-cj 2 m eff a 2 + k 2 a 2 + fc c («2 - a{) + ira eff cJ7 2 a 2 = (S4) 

which can be rewritten as 

-cj 2 m e ff + im e ff^7i \ ^ / fci + k c —k c 



, • f - 1 \ / 7 , 7 , a = (S5) 

Vta + Ka = (S6) 

using a = ( ai ) and by defining the two matrices 11 and K. By solving the eigenwert problem 



a 2 



fi~ 1 Ka = Aa (S7) 



one can obtain the two (rather complicated) analytical solutions. The real part of these solutions contains the 
frequencies while the imaginary part describes the damping constants of the two branches. 



Fitting the measured anticrossing 

These two solutions are then simultaneously fitted to the measured frequencies and Q factors. We approximate 
k n = fco + K n (U — Uq) + X n (U — /To) 2 . This is a second order series expansion of the parabolic frequency tuning [I4j[26] 
around the crossing voltage Uo- The fit shown in Fig. 3 yields the following values: fco = 3.2 ^, k\ = 9.3- 10 -3 k 2 = 
-12.5 • 10- 3 Ai = 0.41 • 10- 3 A 2 = -0.57 • 10" 3 fe c = 3.76 • 10" 3 ^ and U = 6.547 V using an effective 
mass of m e ff = 1.85 • 10 _15 g. As \U — Uq\ is always less than 3 V in the experiment, the maximal relative influence of 
the quadratic term \ u ~ u °\ Xn [ s below 15 %. 

The fit in Fig. 3 shows only the measured frequencies, as the quality factors of the two modes are nearly identical. 
A resonator on a different chip which was tuned using the microwave power instead of the DC voltage [26 , exhibited 
higher, dissimmilar, quality factors and thus allowed to simultaneously fit frequencies and quality factors. The 
resulting graphs are shown in Fig. [Si] One can clearly see how the quality factors of the red and blue branch change 
as the system is tuned through the coupling region and the oscillation transforms from an in plane to an out of plane 
motion and vice versa. 



Analyzing the time- resolved data displayed in Fig. 4 

For each ramp time r, we measure the time-dependent power spectral density at points A and D of Fig. 2 with a 
bandwidth of 1 kHz to have sufficient time resolution. The measurement is started at t = r, but the data used for the 
analysis is the one taken after time 5 = 3 ms to avoid transient spikes in the measurement (as illustrated in Fig. 4a). 
The resulting exponential decays are then fitted using 6' A ' jD (t,r) = 6 , 4 ' jD (r)e _7 ^ + ^Noise? yielding the noise floor 
^Noise, the damping constant 7, identical to the one determined from spectral measurements, and the mode energy 
at time r Sq ,D (t) at points A and D. 
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FIG. SI. Avoided crossing data and corresponding fit of a second resonator tuned via the microwave power, showing the 
coupling behavior of both frequency and quality factor. As the resonator modes are swept through the coupling region, the 
two branches transform from the high quality factor of the out of plane mode to the lower quality factor of the in plane mode 
and vice versa. 



Converting the measured data into transition probabilities 



The transition probabilities of a classical Landau- Zener- Transition [30] are 

^Pdia = e aa and P adia = 1 - P dia . (S8) 
By rewriting the change of the frequency difference between the two unperturbed states 0J n (U) — * I k o-\-k n (U) 



a = d{u,i-W2) = Mgi) - u, 2 {Uj)) - (oMUf) - u^iJJj)) = Aw (g9) 
dt r r 

Pdi&b can be expressed as a function of the ramp time r and the measured frequency differences Auj between the initial 
voltage Ui and the final voltage [//, as the two voltages are kept constant throughout the experiment and only r is 
varied. By introducing the mechanical damping term e~ 7t at t = r, the probability is transformed into a normalized 
state population 

Sdia(r) = e^^ T (S10) 
Sadia(r)= (l-e^)e-^. (Sll) 

The amplitudes of the two datasets Sq(t) and Sq(t) are then rescaled to fit these two equations. Both are shown in 
Fig. 4b of the main text along with the theoretical curves given in (S10). Note that T, 7 and Au are already known 
from previous measurements and are no fit parameters. 



Coupling mechanism 

Hybrid mode shapes 



To learn more about the coupling and to understand how the transformation from an in-plane to an out-of-plane 
motion (or vice versa) takes place during an adiabatic transition, it is interesting to look at the spatial mode profiles 




FIG. S2. The polarizations of the hybrid modes in the coupling region are sketched as black lines in the anticrossing diagram 
shown in (a). A horizontal line represents the in-plane and a vertical line represents the out-of-plane mode. The two hybrid 
modes exhibit polarization directions rotating in the plane perpendicular to the resonator. A mass on a spring model of the 
resonator with an asymmetric "coupling spring" is displayed in (b). 



in the anticrossing region. The solutions of the differential equations 
out of phase (lj = 



SI 



and 



S2 



are an in phase (uj 



k n 
m eff > 



and an 



*m 2 f f c ) com bination of the fundamental mechanical modes. The mode polarizations and their 

showing the transition between the 



a 



qualitative evolution throughout the coupling region is sketched in Fig. S2 
pure in- and out-of-plane modes via the diagonal hybrid modes. 

These two diagonal hybrid modes have different frequencies and thus different energies. In a perfectly symmetric 
beam with a rectangular cross section one would not expect any difference between the two diagonal modes, thus 



the coupling between the modes has to be connected to some asymmetry. This is visualized in Fig. S2^b): the two 
springs labeled k\ provide the restoring force of the out-of-plane mode, the springs ki correspond to the in-plane 
mode. One can directly see from the schematic that the coupling spring k c introduces an asymmetry into the system. 
One coupling mechanism, related to the asymmetric beam position between the electrodes, is discussed in the next 
section. 



Electrical field coupling 




FIG. S3. The red lines of equal electric potential between the two electrodes (calculated using COMSOL finite element 
simulation) demonstrate the inhomogenous electric field (blue arrow) in x as well as z direction. 



A mechanical resonator chip with different gaps between the side electrodes and each beam exhibits different 
coupling constants T for each resonator. The smallest gap (of about 60 to 70 nm) yields T = 7.77 kHz (as presented 
in the main paper). Fitting the frequencies of the other beams in their respective coupling region gives a T of 
7.27 ± 0.09 kHz, 6.10 =b 0.02 kHz and 5.31 ± 0.03 kHz with increasing gap size up to roughly 150 nm. As the electric 
field between the electrodes decreases with their increasing separation d (and the voltage is approximately constant), 
the coupling seems to be mediated by the electric field, even though the exact gap sizes are unknown and thus a 
quantitative relation can not be established. 

This electrical coupling between the in-plane mode (oscillating in x direction) and the out-of-plane mode (oscillating 
in z direction) can also be shown using the following simple model: 



Starting from one undamped coupled equation (a simpler version of equation (SI)) 

F x = k x x + k c {x- z), (S12) 
the coupling constant is just the derivative of F x in z-direction: 

a -t - < sis > 

The electric force on the dielectric beam is the gradient of its energy W in an external electric field E 

Fa = -W = -V(p • E) = -V(aE 2 ) (S14) 

using the polarizability a. Thus, the derivative of the x component of F e \ in z direction yields a dielectric coupling 
term 

d_ (_da&\ = _J^_ = _^ ^ 
dz \ dx J dzdx c ' e 

As there is a gradient of the electric field in z direction (the electrodes are above the beam) and in x direction (from 
asymmetry, otherwise the in-plane mode would not tune with the applied DC voltage), fc c?e j is not zero and contributes 



at least partially of the observed coupling strength. This is visualized in Fig. [S3] if the beam is not perfectly aligned 
between the two electrodes, the resulting effective electric field exhibits a gradient in x and z direction. The field- 
dependent coupling mechanism also explains why the data shown in Fig. [Sl| measured with a DC voltage of V, 
exhibits a weaker coupling of less than 2 kHz (the microwave field used to detect the beam motion also leads to an 
electric field, but the effective voltage is smaller). As the two modes can not be tuned into resonance without applying 
a (DC or microwave) electric field, it is not possible to test if there is also any purely mechanical coupling, e. g. caused 
by interactions between the modes mediated by the prestress of the beam or coupling effects in the shared clamping 
points of the two modes. 



